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Abstract 

A  time  dependent  definition  of  the  virtual  height  of  a  reflected  wave 
train  is  suggested.  Ihis  definition  is  such  that  its  accuracy  increases  as 
the  width  of  the  incident  pulse  increases.  Moreover  it  is  theoretically- 
possible  to  obtain  an  estimate  of  the  virtual  height  with  as  small  an  error 
as  desired  no  matter  what  the  nature  of  the  reflecting  medium. 

Suppose  an  incident  pulse  having  width  W  and  a  carrier  frequency  co^ 
produces  a  reflected  wave  which,  when  measured  at  a  fixed  point  x  =  0  in  space, 
is  R(co  ,W,t).  The  suggested  definition  of  the  virtual  height  h'(co^)  is 


W{ai^)   =  lim  c  • 
w->oo 


00  /   00 


(      tlR(co^,w,t)l^  dt/  f  lR(co^,w,t)|^dt 

-00  -00 


00  /  °°  '\ 

-  [  t|l(co^,w,t)|^  dy     [  Il(co^,w,t)l2dt  I  , 

-00  -00  J 


where  c  is  the  free  space  velocity  of  light.  This  relation  for  fa'C^i^Q)  holds 
for  any  physically  reasonable  incident  wave  train  and  is  independent  of  its 
envelope  shape. 

An  example  is  given  of  a  reflected  wave  whose  virtual  height  cannot  be 
determined  by  inspection  in  the  usual  manner.  The  expression  given  here  for 
the  virtual  height  provides  the  correct  result. 
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1.  Introduction 

Much  of  the  data  obtained  from  radio  sounding  investigations  of  the  upper 
atmosphere  is  given  in  terms  of  the  virtual  height  reached  by  a  narrow  pulse- 
modulated  wave  transmitted  from  the  earth.  A  virtual  height  record,  which  is  a 
record  of  the  group  delay  time  of  the  test  pulse  multiplied  by  the  free  space 
velocity  of  light,  is  generally  made  from  measurements  of  the  time  required 
roughly  for  the  wave  packet  to  go  out  from  and  return  to  the  antenna  from  which 
it  was  radiated. 

Some  ambiguity  in  the  measuranent  of  the  delay  time  of  the  pulse  will  always 
occur  in  a  standard  experiment.  The  reason  for  this  ambiguity  lies  in  the  difficulty 
of  prescribing  a  fixed  point  of  reference  on  the  pulse  envelope  so  that  the  move- 
ment of  the  wave  group  can  be  described  as  the  movement  of  a  single  entity.  The 
difficulty,  of  course,  occurs  because  the  pulse  envelope  must  undergo  a  certain 
amount  of  distortion  in  its  travel  back  and  forth,  including  some  spreading  of  its 
energy  over  a  longer  time  interval  than  that  which  determined  its  inital  width. 

Still,  in  many  cases  the  amount  of  distortion  suffered  by  the  pulse  is  small ^  J . 
When  such  a  case  is  under  consideration  the  ambiguity  in  the  virtual  height  measure- 
ment should  be  coirespondingly  small.  It  may  be  useful,  however,  to  remark  that 
in  any  case  the  ambiguity  caused  by  the  distortion  of  a  test  pulse  in  a  radio  sounding 
experiment  is  not  inherent  in  this  kind  of  measurement.  It  is  due,  rather,  to  the 
imprecise  concept  of  group  delay  time  used  in  the  interpretation  of  the  data  supplied 
by  the  experiment.  In  the  present  paper  an  attanpt  will  be  made  to  overcome  the 
difficulty  by  referring  to  the  precise  definition  of  group  delay  time  instead  of  the 
heuristic  idea  of  it. 

The  precise  definition  of  virtual  height  is  generally  given  as  the  free  space 
velocity  of  light  times  the  derivative  with  respect  to  frequency  of  the  phase  of 
the  complex  reflection  coefficient  which  is  characteristic  of  the  medium  being 
investigated^  J » L  J,  To  be  more  explicit:  imagine  a  c.w.  of  unit  amplitude  and 
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frequency  co  incident  on  the  mediuir.  The  complex  number   |r(oi)|   expn^(a))j giving 
the  amplitude   |r(co)|   and  the  phase  j((co)  of  the  reflected  wave^is  the  reflection 
coefficient  of  the  medium;  the  virtual  height  of  a  pulse  at  a  frequency  oo    for  this 
medium  will  be 

(1.1)  h'   -  c  d^/dojj 


CO    »    CO 
0 


where  c  is  the  free  space  velocity  of  light. 

If  this  precise  definition  of  virtual  height  is  used,   a  second  criticism  of  the 
standard  measurement  technique  comes  to  mind.     Due  to  the  so-called  uncertainty 
principle  relating  a  transient  wave  form  and  its  frequency  spectrumja  very  narrow  pulse 
will  have  a  spectrum  which  is  rather  evenly  distributed  over  a  broad  frequency  band; 
consequently,  the  idea  of  the  group  delay  time  of  the  pulse  loses  its  intuitive  meaning. 
The  group  delay  time,   then,  will  not  necessarily  be  given  correctly  in  terms  of  the 
phase  derivative  of  the  reflection  coefficient,  at  least  according  to  the  arguments 
ordinarily  found  in  literature*- J *'--'*  L^-^.     In  fact,   in  order  to  obtain  a  sensible 
definition  of  group  delay  time,   one  must  deal  with  a  wave  packet  having  a  narrow-band 
frequency  spectinim,   that  is,   a  wave  packet  considerably  spread  out  in  time.     Therefore, 
it  seems  reasonable  that  in  order  to  improve  the  accxiracy  and  to  remove  the  ambiguity 
of  virtual  height  measurements,   one  should  use  as  broad  a  pulse  as  possible  in  radio 
sounding  experiments. 

It  is  the  purpose  of  this  paper  to  suggest  the  kind  of  measurement  appropriate 
to  the  use  of  a  broad  test  pulse  instead  of  a  nairrow  one.     The  analysis  shows  that 
theoretically,   at  least,  the  ambiguity  in  such  a  measurement  can  be  made  arbitrarily 
small    for  all  carrier  frequencies  no  matter  what  the  dispersive  properties  of  the 
reflecting  medium  may  be.     In  particular,  the  h    -  f  curves  obtained  at  the  present 
time  from  radio  soundings  of  the  ionosphere  exhibit  a  great  deal  of  ambiguity  when 
the  carrier  frequency  is  near  the  penetration  or  critical  frequency  of  a  layer.     Thus 
the  use  of  a  broad  test  pulse  with  the  interpretation  of  virtual  height  suggested  here 
offers  the  possibility  for  improving  the  accuracy  of  measurements  at  frequencies  in 
the  neighborhood  of  the  critical  frequency,  for  example. 
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2*     Definitions 

The  proposed  modification  of  the  -method  for  obtaining  the  virtual  height  is  based 

on  the  following  considerations. 

Suppose,  at  a  point  x  ■  o  in  free  space  below  the  medium  to  be  investigated, 
the  incident  pulse  is  given  by  a  real  function  of  time  t 

(2.1)  I(cOjj,  w,t)  •  I(t/w)  cos  0)^  t, 

where  I  (t/w)  represents  a  more  or  less  arbitrary  pulse  envelope  with  a  maximum  at 
t  -  0.  The  constant  w,  assumed  positive,  roughly  fixes  the  width  of  the  pulse  in 
the  sense  that  as  w  becomeslarge  I (t/w)  remains  near  its  value  at  t  ■  0  for  a  large 
range  of  t.  For  example,  we  might  have 

(2.2)  I  (t/w)  -  exp  (-  tVsw^)  , 

a  Gaussian  envelope  with  standard  deviation  w  or 

f  1    for  -w/2  <  t  <  w/2 

(2.3)  I (t/w)  .  {         ... 

1  0    for  1 1 1  >  w/2 

a  square  pulse  of  width  w. 

If  we  define  the  Fourier  transfonn 

00 

i(co)  -  [    I(t)  exp(ja)t)d^,  l(<o)  -  |i(co)|  exp  |^3©(co)] 
/-oo 

(2.I4) 

I(^)  -  --     i(co)  exp  (-jcot)  dm  , 

/-CD 

then  the  Fourier  transform  of  the  incident  pulse  (2.1)  is  the  function 

00 
^(^%t  ^t^^   ■  I    ^("o*  ***^  ®^  (ja5t)dt 
(2.5)  ^-°° 

-  (W/2]ii{w(cD  +  ©^)|  +  i  (w(co  -  a^)}!      . 
In  terns  of  the  reflection  coefficient  r(<)o)  the  Fourier  transform  of  the  reflected 
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wave  packet  R(i»  »  w,t)  is  r(co)  i(w  »  w,co)  and  we  have 


00 


(2.6)     R("o»  ***)  ■  57     '^("^  i('^  *  w»<»)  ®^  (-jwt)db. 

00 

It  will  be  shown  later  that  under  these  conditions; 

00  /    00 


lin 

W->00 


00  /    00 

(        |R((o^,  w,t)|^  dt/f        I(co^,  w,t)|2  dt  -   Ir(co^)r 

/-CD  /-QD 


•00  '-00 

and 


(2-7)  ,00  y  00 

lim     cjl         t|R(co^,  w,t)|2  dt/j        lR(co^,w,t)  |^  dt 

[/-OO  -^      '-00 

,00  /  00 


1-00 


t|l(cO^,W,t)|^    dt/j  |l(cD^,W,t)|2    dt| 


c  dffM/dcoi^^ .  ^  -  h'  (CO )        ^^y  i-i))- 


O  '   o 

The  relations   (2»7)  indicate  that  except  for  quantities  of  infinitesimal  order  as 
w  becomes  large  the  amplitude  of  the  reflection  coefficient  and  the  virtual  height 
of  the  incident  pulse  can  be  expressed  in  terms  of  the  power  and  the  centroid  of 
the  power  in  the  incident  and  reflected  wave  packets.     The  error  committed  in  using 
(2.7)  with  large  but  finite  w  would  sean  to  depend  on  co  w  as  well  as  on  the  ratio 
of  w  to  the  effective  width  of  the  reflecting  medium. 

In  the  case  of  a  narrow  incident  pulse  which  produces  a  narrow  reflected  pulse 
the  measurement  of  virtual  height  obtained  by  the  standard  technique  (consisting  of 
a  direct  physical  observation  of  the  delay  time)  will  give  the  same  result  as  the 
approximation  based  on  (2.7)  where  the  actual  finite  value  of  w  is  used  instead  of 
w  -«»oo.     Thus,  the  approximation  suggested  by  (2.7)  will  always  be  at  least  as 
accurate  as  the  standard  method  for  measuring  virtual  height  and  should  be  more 
accurate  when  a  broad  incident  pulse  is  used  instead  of  a  narrow  one. 
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3.  Example  of  a  dispersed  reflected  pulse 

It  is  easy  to  find  reflection  coefficients  in  connection  with  a  given  incident 
pulse  of  any  desired  width  such  that  an  inspection  of  the  corresponding  reflected 
wave  packet  will  reveal  no  clue  to  the  virtual  height  as  defined  by  (1«1)*  If  the 
incident  pulse  is  sufficiently  broad,  however,  the  computation  suggested  by  (2«7) 
will  always  give  the  virtueOL  height  even  in  this  case* 

Consider  for  example,  a  reflection  coefficient 

(3.1)  r(6))  -  ja  Qxp  (2J<aL)/(«  ♦  ja)  , 

where  o  and  L  are  real  positive  constants.  The  amplitude  of  r(a))  is 

(3.2)  |r(co)I  -  a/(a)2  *  a^)!/^  , 
the  phase  of  r(co)  is 

(3.3)  ^(co)  ■  -tan"-"-  (o/o))  ♦  n/2  +  2ooL  , 
and  thus  the  virtual  height,  according  to  (1*1)^  is 
(3.U)  h'  (oog)  -  ca/(co^  +  o^)  +  2dL  . 

Let  us  choose  for  the  incident  pulse 
(3.5)  I(a)Q,¥,a))  -  I(t/w)  cos  to^t  , 

where 

{i,      for  Itl  <  w 
0,     for  |t|  >  w 

The  reflected  wave  packet  is  given  by  the  e^qjression 
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10,  t<2L-w 

(3.7)  R(wQ,w,t)  -A  sin  0[exp(a(2L-w^t)]  sin(co^w-©)+3in{cOg(t-2L)+«}]  ,2L-w<rt<2L+w 

sin  9texp(a(2L-»».t)jsin(co  w+9)+exp{a(2L-i»-t)}sin(cOgW-eri,2L+w<t   , 


2         2  1/2 
where       sin  ©  ■  a/(o     +  o)^)  ' 


REFLECTED    PULSE      R(aj^,W,t)    VS     TIME    t 


a=  I 

W=l 
L=5 


o 


The  envelope  of  R(co  ,w,t)  has  essentially  the  same  shape  for  any  value  of 

the  carrier  frequency  co  .     Thus,    there  is  no  obvious  way  to  obtain  the  virtual 

height   (3«U)  from  R(a)  ,w,t)  by  inspection. 

We  can  apply  (2.7)  to  compute  the  virtual  height,   however.     We  havet 

1 
(3.8)  i(o3)  -  I       exp  (Jcot)dt  -  2  simo  /u 


-1 


30  that  by  (2,5) 


(3«9)  i(co  ,w,co)  ■  sin|w(a)+co  )[/(a>+a)  )  +  sinlw(6j-<d^)>/(&j-w  ) 


We  also  havet 


00 


'-00 


|l(co  ,w,t)p  dt  «  w  +  <J(w) 


(3.10) 


.00 


y-00 


|r(co  ,w,t)|      dt 


WiV(a^   *  cob  •>■  (5tw) 


00 

(' 

'-00 


t|R(w  .w,t)|^  dt  «  wtiV(o)^  *  a.^)^  *  2wLaV(a^  +  "^)  +  *"(«)      - 


The  general  results  (2«7)   are  thus  verified  in  this  special  case.     The  error 
actually  depends  on  co  w  and  aw  but  not  particularly  on  L,   as  a  more  detailed 
calculation  will  show.     The  quantity  L,  according  to   (3»7)»  is  just  the  distance 
of  the  edge  of  the  reflecting  medium  divided  by  c,  the  free  space  velocity  of  light, 
The  quantity  a  may  be  regarded  as  a  measure  of  the  width  of  the  reflecting  medium 
since  r(ai)  becomes  zero  if  a  approaches  zero. 


U.     Proof  of  (2.7) » 

Relations  and  definitions  (2»U),  (2.5)  and(2j6)  are  necessary  for  the  proof  of 
(2»7).  A  standard  argument  which  is  employed  in  the  proof  of  the  Fourier  integral 


-  8  - 


theorem  will  be  used  to  demonstrate  a  preliminary  result  we  shall  require. 

Theorem  A«  Let  f  (co)  be  continuous  and  uniformly  bounded  for  all  real 
CO,  and  assume 


/  —00 


ll(oi>)|     dco  <oo     . 


-00 


Then 


00 


'-00 


:{' 


00 

f(a))|i(Ug,w,a))|^  dco  -  Jjf(co^)+f  (-«^)l|        |i((o)|^dio  ♦  «r(w) 


Proof:     By  (2.?)  we  have 


(U.l) 


|i(co^,w,«)|2=.  f^[|i{w(a,fa)^))|^|i{w(a«o^)]  \^*i[Mc^<,^)j±*(v{c^^)j 
*  i*[w(co+cOp)ji^(<o>.ojQ)j       . 


First  we  consider  the  first  two  terms  on  the  right  of   (U«l)  at  the  same  time. 
We  have 


00  00 

I        |f  (a))-f  (ta)^)j|i(w(co+co^)j|2dco-i  j        j^f(67w+(o^)-f( +0)^)1 
/-CD  /-OO 


i(<r)|^d<r 


(U.2) 


rr 

r^ 

+ 

w   J 

-00. 

^  (<5/w+co^  )-f  (*coJK  \±M\  ^do-. 


Consider  the  first  terra  on  the  right  of  (l4.*2)*     We  have 

/S  00 

(U.3)  i  (         |f(o/wico^)  -  f(+<o^)I|i(o-)|V<  i  f       iiCo^I^do- 

for  w  large  enough,  by  the  contin\iity  of  f  (co)  at  co  •  The  last  two  tenns  on  the  right 
of  (U,2)  are  evidently  equal  to  (l/w)  dtl).  Since  e  in  (U»3)  can  be  made  as  small 
as  we  please  for  large  w,  we  have  finally 
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00  00 

(       f(co)|i/w(co;^^)]|^(kB  =    (l/w)f(+co^)|        |i(o-)|^dcr+(l/w)(y(l)      . 
/-OD  /-oo 


•  00  '-00 

Now  we  consider  the  last  two  terms  on  the  right  of   (U.l)  one  at  a  time.     We  have 


.00 


f  (co)i{w(w+a>^)]i   /w(6)-oj^)(dco 


-00 


(U.?)       -  (1/ 


'w)  j     f((r/w+«^)i' 
J-oo 


((T)i({f+2co^w)d<r 


-00 


J 


(1/w)  f((3-/w+co^)i  (a^iCC+aw  w)dcr 


(1/w) 


1  /v         '-00 


f  (<r/w+co^)i   ((5)i(<r>-2w  w)dd^ 

O  0 


By  Schwartz's  inequality  the  first  term  of  the  right  side  of  (U.5)  is  of  order 


[1. 


(1/w)     I         |i(6^rdcri 

'00 


1/2 


-   jZa,  w  +/^ 


.  ■'2co^w  -/w 


1/2 
iKff^l^dd  -   (l/w)<f(l)     . 


The  second  and  third  terms  on  the  right  side  of   (U.5)  can  also  be  estimated  by 
Schwartz's  inequality; 

00 

(1/w)  I       f  (<r/w+a)^  )i*(cr)i  (<5*2co^w)  dsr 


'/w 


(1/w) 


00 


,00 


|i(or)|2dor|      lf(<r/w^^)|2|i(<j-+2<o^w)|2d<jj  '    -(l/w)<y<l)  , 


1/2 


and  of  course  we  have  tJie  same  result  for  the  other  term.     Combining  these  results 
with  (U«l)  we  arrive  at  Theorem  A. 

From  Plancherel's  theorem  for  the  Fourier  transform  (2»6)  and  the  fact  that 
|r(oj)|   is  an  even  function  of  co,  we  obtain  the  first  part  of   (2«7)»     To  obtain  the 
second  part  we  start  with  Plancheir«l ' s  theorem  applied  to  the  first  moment  of  a 
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function.  Thus,  if  £((£>)  ±a   the  Fourier  transform  of  F(t)  we  note  that 


(li.6) 


00 

I   * 

/-oo 


t|F(t)|     dt  • 


00 

(J 


^  I       f  ■    (a))f^(a>)cLo 


Applying  (U«6)  to  the  numerator  of  the  first  term  on  the  left  of  the  second  equation 

in  (2.7)  we  have 

oo  00 

I       t|R(co^,w,t)I^  dt  -  ^  I       ^^r(oo)i((Og,w,to}»  r*(co)i*(«Q,w,a))dcD 
''-00  /-oo 


if         d 

J  "CO 


|!;2f(co)+8(®^,w,a>)}|i(coQ,w,<)o)|^   |r(co)|     dco 


(U.7) 


*  2n  I         3S 
J  -oo 


||r(«)I|i(tD^,w,cc)|j  [^|r(a3)|  |i(co^,w,(o)|Jdffl 


00 
'  -00 


for 


f  ^  [|r(o))|  |i(«Q,w,«)|j  [|r(co)l|i(«.Q,w,ec)|]  do 

/ -00 


/ -00 


A  similar  consideration  applies  to  the  numerator  of  the  second  term  on  the  left 
of  the  second  equation  in  (2«7)«     Combining  these  results  with  an  application 
of  Theorem  A,  we  obtain  the  second  relation  in  (2»7). 
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